hep-ph/970235^ CERN-TH/97-25 

NTZ 97/6 



Scale dependence of the chiral-odd twist-3 
distributions hi{x) and e(x) 

§5 ! A.V. Belitsky 

5— I , 

C3 . Bogoliubov Laboratory of Theoretical Physics 



in 



43 



X 



Joint Institute for Nuclear Research 
14 1980, Dubna, Russia 



> 

in ■ D. Miilleif] 

cs 

o 



Theory Division, CERN, 



CH-1211 Geneva 23, Switzerland 



Abstract 

We evaluate the complete leading-order evolution kernels for the chiral-odd twist-3 distribu- 
tions e(x) and hi(x) of the nucleon. We establish the connection between the evolution equations 
in light-cone position and light-cone fraction representations, which makes a correspondence be- 
tween the non-local string operator product expansion and the QCD-inspired parton model. The 
compact expression obtained for the local anomalous dimension matrix coincides with previous 
calculations. In the multicolour QCD as well as in the large-x limit the twist-3 distributions obey 
simple DGLAP equations. Combining these two limits, we propose improved DGLAP equations 
and compare them numerically with the solutions of the exact evolution equations. 

CERN-TH/97-25 
February 1997 



1 Permanent address: Institut fur Theoretische Physik, Universitat Leipzig, 04109 Leipzig, Germany 



1 Introduction 

The factorization theorems |J provide us with powerful tools to study the processes with large 
momentum transfers. They give us the possibility to separate the contributions responsible for 
physics of large and small distances. The former are parametrized by the hadron's parton distri- 
bution functions, or by parton correlators in general, which are uncalculable at the moment from 
the first principles of the theory, while the second ones — hard-scattering subprocesses — can be 
dealt perturbatively The parton distributions are defined in QCD by the target matrix elements 
of the light-cone correlators of field operators 0. This representation allows for the estimation 
of these quantities by the non-perturbative methods presently available, which are close to the 
fundamental QCD Lagrangian ||. 

The increasing accuracy of the experimental measurements requires the unravelling of the twist- 
3 effects in the hard processes, which manifest the quantum mechanical interference of partons in 
the interacting hadrons. The most important advantage of the twist-3 structure functions is that, 
while being important for understanding the long-range quark-gluon dynamics, they contribute 
at leading order in 1/Q (Q being the momentum of the probe) to certain asymmetries 0, [J 
and therefore may be directly extracted from experiments ||. To confront the theory with high- 
precision data, the knowledge of the scale dependence of the measurable quantities is needed. 
The Q 2 evolution of the parton distributions [0, |^| can be predicted by exploiting the powerful 
methods of renormalization group (RG) and QCD perturbation theory. The evolution equation 
of the twist-3 polarized chiral-even nucleon structure function gi{x) was extensively discussed in 
the literature [p|— 1|13| as well as its solution in the multicolour QCD and asymptotical values of 
orbital momentum n — > oo (x —>■ 1 region) |H| . 

In this paper we address the study of the evolution of other twist-3 structure functions of the 
nucleon: chiral-odd distributions e(x) and hi(x) ||, which open a new window to explore the 
nucleon content. The present study was induced by several reasons and pursues various goals. 
First, while the anomalous dimensions of the local operators corresponding to the moments of the 
chiral-odd distributions were computed for the real QCD case |T5|] , the kernels of the evolution 
equations were derived only in the multicolour limit |16[. From the point of view of experimental 



measurements, the knowledge of the evolution of the whole x-dependent distribution is welcome. 
However, exact equations that takes into account the 0(1/N 2 ) effects are lost, although we may 
expect sizeable 1/N 2 effects for the small- x behaviour of the distribution functions, provided the 
non-leading terms in N c yield the rightmost singularity in the complex n-plane of the angular 
momentum with respect to the leading-order result. Second, up to now the relation between 
different formulations of the evolution equations in the light-cone fraction ||, 10] and light-cone 



position [IT], [L7| representations was obscure. So, the aim of the present paper is two fold: to fill the 



gap in our knowledge of the exact (taking into account of the 1/N 2 effects) higher-twist evolution 
equations for chiral-odd distributions and to provide the relation between different calculational 
approaches. We attempt here to clarify these issues. 

The outline of the paper is the following. We start with the construction of the basis of the 
twist-3 chiral-odd polarized and unpolarized correlation functions, mixed under the renormaliza- 
tion group evolution, and find the equations that govern their scale dependence in the momentum 
fraction as well as in the light-cone position representations in Section 2. In Section 3 we present 
the Fourier transformation of the evolution kernels, which simplifies the transition from one repre- 
sentation to another. Section 4 is devoted to the construction of the generalized DGLAP equation 
in mixed representation for the twist-3 distributions, which is the most useful for solving the RG 
equations in the large- N c limit as well as for asymptotical values of the orbital momentum n. The 
corresponding analytical solution, as well as a numerical study of exact equations is performed in 
Section 5. The final Section contains a discussion and concluding remarks. 

To make the discussion complete and more transparent, we include several appendices. In 
Appendix A we present definitions and some properties of the B-functions that appeared in the 
formulation of the evolution equations in the momentum space. In Appendix B the evolution 
equations for the redundant basis of operators are found in the Abelian gauge theory. This shows 
the self-consistency of the whole approach. The anomalous-dimension matrix of local twist-3 
operators obtained from the evolution equations we have derived in the main text are written in 
Appendix C. It coincides with results known in the literature p 



2 Evolution of chiral-odd twist-3 correlation functions 

As we mentioned in the Introduction, the parton distribution functions in QCD are defined by the 
Fourier transforms along the null-plane of the forward matrix element of the parton field operators 
product separated by an interval A on the light cone: 

:F(A)=.F(A,O) = 0*(O)0(An), (1) 

where <fi denotes a quark j/iora gluon field B^. We suppress the dependence on the renormalization 
scale /j,r, necessary to make this equation well defined in field theory. Throughout the paper we 
use the ghost-free B + = B^n^ = gauge. Here n is a light-cone vector n 2 = normalized with 
respect to the four-vector P = p + \M^n of the parent hadron h of mass M^, i.e. nP = 1, and 
p is a null vector along the opposite tangent to the light cone such that p 2 = 0, np = 0. In 
any covariant gauge a path-ordered link factor should be inserted between the 0-fields so as to 
maintain the gauge invariance of the physical quantities. The Fourier transformations from the 



coordinate to the momentum space and vice versa are given by 

F(x) = [p e iXx (h\F(\)\h), (h\F(\)\h) = fdxe- iXx F(x). (2) 

Both of these representations display the complementary aspects of the factorization. The light- 
cone position representation is suitable to make contact with the operator product expansion 
(OPE) approach, while the light-cone fraction representation is appropriate for establishing the 
language of the parton model. Throughout the paper we will use the light-cone position and the 
light-cone fraction representations in parallel. 

The multiparton distributions corresponding to the interference of higher Fock components 
in the hadron wave functions that emerge at the twist-3 level are the generalizations of flip to 
3-parton fields 

T{\, fi) = F{\, 0, fi) = 0*(/in)0(O)0(An). (3) 

We do not display the quantum numbers of the field operators since they are not of relevance at 
the moment. The direct and inverse Fourier transforms are 

F(x,x') = [—^e iXx -^ x '(h\F(\,i2)\h}, (h\F(\,fi)\h)= fdxdx'e- iXx+i ^'F(x,x'). (4) 

The variables x and x' are the momentum fractions of incoming <ft and outgoing cff partons, 
respectively. The restrictions on their physically allowed values come from the support properties 
of the multiparton distribution functions discussed at length in Ref. [19|], namely F(x, x') vanishes 
unless < x < 1, < x' < 1. 

Beyond the leading-twist level the intuitive parton-like picture is not so immediate, as one 
usually starts with an overcomplete set of correlation functions. However, the point is that the 
equations of motion for field operators imply several relations between correlators, and the prob- 
lem of construction of the simpler operator basis is reduced to an appropriate exploitation of these 
equalities. The guiding line to disentangle the twist structure is clearly seen in the light-cone for- 
malism of Kogut and Soper ]n|. Consider, for instance, the correlators containing two quarks ipip. 
Then decomposing the Dirac field into "good" and "bad" components with Hermitian projection 
operators V± = |7^7±: ip± = y±ip, we have three possible combinations ip\.ip+, ip+ip- ± w-?P+, 
and ip_ip-, which are of twist 2, 3 and 4, respectively. The origin of this counting lies in the 
dynamical dependence of the "bad" components of the Dirac fermions 

tp- = - l -d- 1 (tP ± + m) 1+ tlj + . (5) 

These components depend on the underlying QCD dynamics, i.e. they implicitly involve extra 
partons and thus correspond to the generalized off-shell partons, which carry the transverse mo- 
mentum. For this reason we come back to the on-shell massless collinear partons of the naive 



parton model, but supplemented with multiparton correlations through the constraint (H). The 



operators constructed from the "good" components only were named quasi-partonic [p!0| . The 
advantage of handling them is that they endow the theory with a parton-like interpretation for 
higher twists. 

At the twist-3 level the nucleon has two chiral-odd distributions e and h^ Q, which can be 
measured in the polarized Drell-Yan and semi-inclusive DIS processes. Under QCD evolution, 
they couple to complicated quark-gluon operators and correlation functions, depending on the 
quark mass and intrinsic transverse momentum. We will treat below the unpolarized and polarized 
cases separately. 

2.1 Unpolarized distributions 

In the unpolarized case we define the following redundant basis chiral-odd twist-3 correlation 
functions: 

e(x) = y^(hmO)^n)\h), (6) 

M(x) = \ [Qe**»{h\$(0)rnry + 1>{Xn)\h), (7) 

L J ATI 

Di(x,x') = \ f^e^-^'(h\^n)g 1+ ^(0)i;(\n)\h), (8) 

Z J 2,71 Z7T 

D 2 (x',x) = y^ e ^'- lXx (h\i;(\n)g^ ± (0) 1+ ^n)\h). (9) 

The functions D\ and D 2 are related by complex conjugation [Di(x, x')]* = D 2 (x',x). The quan- 
tities determined by these equations form a closed set under renormalization; however, they are 
not independent, since there is a relation between them due to the equation of motion for the 
Heisenberg fermion field operator: 

e(x)-M(x) - fdx'D(x,x') = 0, (10) 

where we have introduced the convention 

D(x,x') = -[D 1 (x,x') + D 2 (x',x)}. (11) 

This function is real-valued and antisymmetric with respect to the exchange of its arguments: 

[D(x, x'))* = D(x, x'), D(x, x') = -D(x', x). (12) 

In Appendix B we present a set of RG equations for the correlation functions determined by 

Eqs. ©-© derived in Abelian gauge theory. The relation (|ToD provides a strong check of our 

2 In the language of OPE the local twist-3 as well as twist-2 operators contribute to the matrix elements of the 
distribution function Kl- Equation ( p7|) is a consequence of this fact. 



calculations^. It allows the reduction of the RG analysis to the study of scale dependence of the 
three-parton D and mass- dependent M correlators only. 

The function D(x,x') is gauge- variant provided we use the gauge other than light-cone; there- 
fore, we are forced to introduce the gauge-invariant quantity 

Z(x, x') — (x — x')D(x, x'). (13) 

Using the advantages of the light-cone gauge, where the gluon field is expressed in terms of the 
field strength tensor (the residual gauge degrees of freedom are fixed by imposing an antisymmetric 
boundary conditions on the field, which allows unique inversion): 

B^Xn) = d- l G +tx (\n) = - / dze(\ - z)G+^z), (14) 

and taking into account the relation 

- / — e ±a *e(A -z) = ±—VY-e ±izx , (15) 

2 7 2vr v ; 2vr x y ' 

we can easily obtain from Eqs. @ and (0) the definition of the gauge-invariant quantities in terms 
of three-particle string operators. Generically 

Z(x,x') = U^e^-^'(h\Z(X^) + Z(-^-X)\h), (16) 

where 

Z(X,fj) = Z(\,0,ii) = ^(fin) g G +p (0)cr^(Xn). (17) 

In the same way, for a mass- dependent non-local string operator 

777 _ 

7vF(A) = M>{\ 0) = -^(0) 7+ (^+(A))^(An), (18) 

the Fourier transform is 

r fIX 

M j (x) = x j M(x) = / — e lXx (h\M j {X)\h). (19) 

For the spin- dependent scattering discussed below, the only difference is that one should insert 

also a 75-matrix between the fields in the definitions of the string operators (f[7|) , (|i~8|) . 

3 This fact follows from general renormalization properties of gauge-invariant operators as one expects that the 
counter term for the equation of motion operator can be given only by the operator itself. Its matrix element, 
being taken with respect to the physical state, decouples completely from the renormalization group evolution. 



2.2 Polarized distributions 



Analogously, the set of correlation functions for the polarized case is as follows: 



h(x) 



-9 X 
2 b ° 



dX , 



i \x 



2ir 



(h\i>(0)ia + ^ 5 i/j(\n)\h), 



h L (x) = | J ^e iXx (h\mi<r + -ls^(\n)\h), 



M{x) 



dX , 



i\x 



2vr 



{h\ip{0)mj + ^ 5 ip{Xn)\h) , 



K(x) = \ [^e^(h\mil+P±lsi>(\n)\h), 



1 f d\ dp &, -,,,,•, 

2 7 27T27T 



Di(x,x') = ^ / ^^e^^(h\$(»n)g>y + x (O)>y B iJ>(\n)\h), 



D2(x',x) = 1 -[^%^'~^ 

I J Z7T Z7T 



(/#(An)g7 + ^(0) T5 V>(/m)|/i> 



(20) 
(21) 
(22) 
(23) 
(24) 
(25) 



where S^r denotes the transverse polarization vector of the hadron h (S 2 = —M%). The derivative 
in the correlation function K(x) acts on the quark field before setting its argument on the light 
cone. 

Besides the identity arising from the equation of motion 

(26) 



(27) 



h L (x) - M(x) - K(x) - / dx'D(x, x') = 0, 

there is an equation provided by Lorentz invariance 

d „, s „ f , ,D(x,x') 



a f L>{: 

2xhi(x) = 2hAx) — x—-K(x) — 2x dx - — 

ox J ix' 



— x) 



It means that both parts of this equality are expressed in terms of matrix elements of different 
components of one and the same twist-2 tensor operator. Again we have introduced the C-even 
quantity D, which has the properties 

D(x,x')]* = D{x,x'), D(x,x') = D{x',x). (28) 

Combining the two Eqs. fl2"6| ) and (^7]) we can obtain the following relation between the correlators: 

Z{x,x'), (29) 



2-x— I h L (x) =2xh 1 (x) -x—M(x) + f dx'—^- 



dx 



dx 



x — x 



8_ _d_ 
dx dx' 



where Z(x, x') is a gauge-invariant quantity introduced by Eq. ([13]), but for the function D(x, x'). 
Solving the differential equation with respect to hi{x) the integration constant can be found from 
the support properties of the distribution: hi(x) — for \x\ > 1. The solution is 



dp 



h L (x) = 2x 2 -£h 1 (p) + M(x)-2x 2 -^M(P) + x 2 -£ 



P 



dp 



P 3 



dp r dp' 



x p 2 J p'-p 



d_ _d_ 
dp'dp' 



z(p,p'). 

(30) 



A similar relation was found by Jaffe and Ji Q in Ref . || . Here the dynamical twist-3 contribution 
is explicitly related to the particular integral of the three-parton correlation function Z . In terms 
of local operators it looks like 



{n + 3)[hi 



2[hi] n+ i + (n + l)M n + £(n - I + l)Z l n , 



1=1 



and the definition of moments of distribution functions is given by Eq. ( |C1| ) . 
As before, excluding the functions fl2"T| ) and (P3"D, and using the relations (| 
chose the basis of independent functions in the form: hi(x), M(x), D(x, x' 



and 



(31) 



, we can 



2.3 Evolution equations 

Note that in the leading logarithmic approximation the evolution equations that govern the Q 2 - 
dependence of the three-particle correlation functions are the same, discarding the mixing with 
the quark mass operator. Therefore we omit the "tilde" sign in what follows. We evaluate the 
evolution equations using the different approaches described in || [K]] and |H| and keep to general 
remarks. 

The peculiar feature of the light-like gauge is the presence of the spurious IR pole l/k+ in the 
density matrix of the gluon propagator: 



D^(k) 



fc 2 + i0' 



d 



fJtU 



9p.v 



k. 



(32) 



The central question is how to handle this unphysical pole when k+ = 0. In the calculation of the 
relevant Feynman diagrams shown in Fig. [I] (and trivial self-energy insertions into external legs) 
we assume two different approaches, which employ the principal value (PV) and Mandelstam- 
Leibbrandt prescriptions (ML) [pCR : 



PV ^7 

ML— 



l 

2 



(fen) + iO 
(kn*) 



(kn) — iO 



(33) 
(34) 



(kn)(kn*) + i& 

with the arbitrary four-vector n* satisfying n* 2 = 0, nn* = 1 (without loss of generality, we can 
put it equal to p). The first prescription will be used in the momentum space P, [TO, plf, the 



second one in the coordinate space formulation [|13j . As a by-product we verify that both of them 
do lead to the same result. 



'Corresponding expressions in Ref. pi contain misprints. 



In the light-cone fraction representation we get for the correlation function D(x,x') 



D(x,x') = - — l-C_ 



a \ „ [x x) [j A/f//-- ^ Q o / / / ^ j.„i//'„aqO 






x'M(x)G° u (x', x' - x) ± xM(x')G° n {x, x - x') 



+[D(J3 + x', x') - D(x, x')} . ^ X ~ X '\, &n(x -x',x-x'- (3) 



(x, x — j3) 



(x — x' — (3) 



(d((3, x')—j C —-Q° ll (x, x - (3) + D(f3 + x', x')e° u (x - x', x - x' - p)\\ 



(P + x- x') 
x' 

+ <C F - ?f) (d(/3, x>) ^ + * ~ X >] e° m {x, x-x',x-x' + P) 

+[D{p,x' - x + p) - D{x,x')}^-Q n {x,x - p) n 

+ J dp' (c F D(x, P')^e° n (x', x>-x) + C f ([D(x, P') - D(x, x')] X '_ Q° u (x', x> - P') 

+ [D(x, P' + x)- D(x, x')] V - x \ Q° u (x' -x,x'-x- P') 

(x' — x — P') 

+ (P' + x'-x) f , x' e o (y / ;X / _ ^ + D( ^/ + a; )e° 1 (s / - x ,x'-x- P') 

x \ (x — x') 



C A 



+ [ C F -=f)[D( X ,fi) 



\~) 1 1 1 1 X ■ Ju Ju « Ju Ju ~\ Lj J 



x 



+[D{x _ x ' + ^ p>) _ D{X} a ^]_^_ e o L ( a / i x > _ pr t 

and for the mass-dependent correlation function we have 



a. 



M(x) = -C F — J dpM{p) { 2 



' Q° u (x,x-P) 



(x-P) 



+ 



, -C F D(x,x f )\, 



P±Ze° n (x,x-p) 



(35) 



(36) 



where we have used the dot as short-hand for the logarithmic derivative with respect to the 
renormalization scale ' = /ifj9/cty/|. and the standard plus-prescription fulfilling J dx[...] + = [for 
a definition see Eq. (|B.4|)]. The kernel of the last equation resembles the non-singlet splitting 



function of (un)polarized scattering up to the additional term ^-^C F S(x — P), which is the one- 
loop renormalization constant of the quark mass taken with a minus sign. An explicit form of 
the 0-functions is given in Appendix A. Throughout the paper the plus and minus signs in the 
mass-operator term correspond to the functions D (for e) and D (for Hl), respectively. 

For the string operators (or their matrix elements) we obtain the following compact RG equa- 
tion: 

Z{\n) = ^J o 1 dyJ o y dzi [ C F f5(z)[M 1 (\-w)±M 1 (\y-Li) 



Ca 
2 



2z+[N(y,z)] + --5(y)8(z) 



(C F -^ 



[L(y,z)] + --S(y)8(z 



2z [Z(—Xy, jj, — \z) + Z[X — fiz, —fiy)} 



Z(Xy, fi- Xz) + Z(X - fiz, ny)} 
Z(Xz + (j,z, \iy + \y) 



(37) 



with 



1 — ' z 

[N(y, z)]+ = N(y, z) - 5(y)5(z) f dy' f dz'N(y', z'), N(y, z) = S(y - z) V - + 8{z) V ~, 

Jo Jo y y 

[L(y, z)] + = L(y, z) - 8{y)S{z) C dy' F dz'L(y', z'), L(y, z) = 8(y)- + 5{z) V -. 

Jo Jo z y 

The equations written so far should be supplemented by the following 

"2 



M\X) 



2^ CF Io dy 



a 



-C F / dy 

Z7T JO 



J + 



-2-y-y 2 \M\Xy), 
-2 + ^5(y)\h 1 (Xy). 



(38) 

(39) 
(40) 



The last one, when transformed to the momentum space using the formulae of the next section, 
coincides with the result obtained in Ref. IF22 . 



3 From light-cone position to light-cone fraction 
representations 

Having at hand the evolution equations in different representations for the same quantities, it is 
instructive to relate the kernels in both cases. Such a bridge can be easily established using the 
Fourier transformation for the parton distribution functions given by Eqs. (@) and (§). 

First, we come to the simpler case of the two-particle correlation functions T . The evolution 
equation in the light-cone position space is of the following generic form 



T{X) = / dyK(y)F(Xy), 



(41) 

where JC(y) is an evolution kernel in the coordinate space. By exploiting the definitions (0) we 
can recast the Fourier transform on the language of two-particle evolution kernels. In this way we 
find the direct transformation 

K(x,P)= j 1 dylC{y)8{x-y(5). (42) 

Jo 



And using the general formula 

J o 1 dyf(y)5(x-yP) = f (^\ &° n (x,x-(3). 



(43) 



with the familiar B-function of the momentum space formulation, we can observe that the RG 
equations for two-parton correlators derived in the previous section indeed coincide. The inverse 
transformation can be done 



^e~^ + ^K(x,(3) = fdylC{y)5^-yX) 
2ir Jo 



with the help of the formula 



dj ^-f (fj ©ii(*, x - 0)e-»*«* = £ dyf(y)5(v - yX). 



(44) 



(45) 



The corresponding transformation for three-particle correlators is a little bit more involved. 
The general form of the evolution equation for the light-cone string operator Z(X,fi) reads 



/•i ry 
Z(X,fi)= dy dzlC(y,z)Z(r] n X + r] 12 fi,r] 21 X + r] 22 fi), 
Jo Jo 



(46) 



where rjij are linear functions of the variables y, z. In the momentum fraction representation the 
evolution equation looks like 



Z{x,x') = / dpdp'K{x,x',P,P')Z{P,p'). 



(47) 

Specifying the particular form of the functions rjij, we list below the corresponding conversion 
formulae. 

For the C^/2 part of the evolution equation, the Fourier transformation gives 



K(x, x', (3, (3') = 8{(3' - x') I dy T dzK(y, z)S(x - x'z - (3y). 

Jo Jo 



(4S 



(49) 



The particular contributions are 

K 1 (y,z) = 6{y)6(z) % K 1 (x,x',P,f3') = S(f3-x)S(^-x > ), 

9 

K, 2 {y,z) = 8{z) V ~ + 8{y-z) V - ^ K 2 (x,x', (3,(3') = 8{(3' - x'){-^e° n (x,x - (3) 
y y 1(3 -x 

^-^^) 9!i "-^'-4 (50) 

JCs(y, z) — 1 — > K 3 (x, x', (3, (3') — 8{(3' — x')Ei(x, x — x', x — (3) 

= -8(P'-x')G° ni (x,x-x',x-(3), (51) 

fCJy, z) = z — > KAx, x',(3,(3') = 8((3' — x')\ E%(x,x — x',x — (3) 

I x' 

+ Ez 2 (x,x-x',x-(3)\. (52) 

x' J 

Here we have used ( |43"D and the following useful formula 

E n (x, x - x', x - (3) = f 1 dyy n e° n ((x -(3) + y/3, (x - f3) - y(x> - (3)) 

Jo 



1 
n 



1 - 



' (3-x" 



e? 1 (x,x-a;') + -4 
nx' 



' (3 — x s 



'(3-x 



Q° u (x,x-P). (53) 



10 



The analogous results for the (Cp — Ca/2) part of the kernel read 



K(x, x', (3, (3') = 8{(3' - x') f dy T dzJCty, z)8(x - x'z + (3y), (54) 

JO Jo 



and in particular 



ICi(y,z)=z — > K\{x,x', (3,(3')= 5((3' — x')< Ei(x' — x, — x, x' — x — (3) 

H — S 2 (x' — x, — x, x' — x — (3) \. (55) 

x J 



In addition, we have 



K(x, x', (3, (3') = [ dy f dzKiy, z)8(x - (3z + f3'y)8(x' - (3'y + (3z) (56) 

Jo Jo 



and 



K 1 {y,z) = 5(z) V - % K 1 (x,x',P,P')=8(P-(P'-x' + x))-^Q° n (x',x'-P'). (57) 
y p' — x' 

These formulae complete the list of transformations. Collecting particular contributions, we can 
easily verify that the evolution equations given by Eqs. ( |35|) and ( p7|) agree with each other. It 
should be noted that it is sufficient to have at hand Eqs. ( fl"3D and (^) to perform the conversion 
from one representation to another. 



4 Generalized DGLAP-type equations 

To study the large-iV c limit and establish the relation to the evolution equation for the par- 
tially Mellin-transformed operators introduced in Refs. |11], [H], [L6j , we proceed to the generalized 
DGLAP representation of the evolution equation for three-particles distributions first given in the 



second paper of Ref. fL3 |. For this purpose we define a new function, Fourier-transformed with 
respect to the A variable only: 

z(x,u) = - f—e iXx (h\Z(u\,-u\)±(u^u)\h), (58) 

2 J 2ir 

which is even under charge conjugation and depends on the variables x and u. The latter has the 
meaning of the relative position of the gluon field on the light cone. For < u < 1 the gluon field 
lies between the two quark fields. Because of the support property \x\ < max(l, \2u — 1|), the 
variable x is then restricted to \x\ < 1 and can be interpreted as an effective momentum fraction. 
The evolution equation for z(x, u) can be derived in a straightforward way from the RG 
equation ([37]) for the non-local string operator Z. It can be presented in the form of a generalized 
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DGLAP-type equation: 
z(x,u) = 



ocs_ r dy 
2tt J y 



mi 



dv<P zz (y,u,v)z [-,v\ + P zm (y,u,v)m (-) 



(59) 
(60) 



Here, m(x) = xM(x) and the integration region is determined both by the support of z(x, u) and 
by the kernels 



P zz (x,u,v) 



C F - 



C t 



0i(x, u, v) [L(x, u, v)]+ — 6*2(3:, u, v)M(x, u, v) — -S(u — v)S(x) 



+ ^-®3(x,U,v) 



7 
M(x, u, v) + [N(x, u, v)\ + — -S(u — v)S(x) 



+ 



4 

'u — > u" 

V ->• V , 



X 



P zm (x,u,v) = Cfx 9(x)9(x) — [S(v — u — xu) ± 5(v — u — xu)] , 



■Lm.m. \ -b I ^- / f 



2 — x — x^ 



(61) 



(62) 



(63) 



where the auxiliary functions are defined by: 



Oi(x, u,v ) 
2 (x,u,v) 

3 (x,u,v) 

L(x,u,v) 
M(x,u,v) 
N(x,u,v) 



8(x)8(u — xv)0(u — xv), 

xv \ „ / 1 — xv \ „ / x — u 
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11 



11 



(*\ 9 (™] 9 (?Lp 



uj \ u 

2 

S(u — xv), 



II 



v(v — u) 
2x(l — xv) 



u 3 
ve(u) 



— o(x) H out — xt> J 

■u 1> 



(64) 



The plus-prescription for the arbitrary function A is defined by the equation 

Oi(x, u, v) [A(x, u, v)] + 

= ©i(x,u,v)A(x,u,v) — 8(x)S(u — v) / dx' / dv' 0i(x',u,v')A(x',u, v'). (65) 

Note that, due to the evolution, the variable u is no longer restricted to the region < u < 1. 
Going further, we introduce the Mellin transforms 



z (u 



dxx n z(x, u) and ml' 



dxx n m(x) = / dxx n M(x), 



(66) 
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where n is the complex angular momentum. Operators with different n do not mix with each 
other and satisfy the evolution equations 



z n (u) 



m 



g J dv {/»(«, v)z"(v) + 5(u - v) R>) ± PIM] >»"} 



gC„m" 



The kernels are given by 



P$ z {u,v)=[C F - 



C A 



©i(«, v) [L n (u, v)]+ - 2 (u } v)M?(u, v) - -8{u - v) 



+^0 3 (u,v) 



M2(u,v) + [N n (u,v)] + -~S(u-v) 



+ 



,t) -» W 



(67) 
(68) 



PlJv) = Q 



2 - v n [2 + n{2 + (n + l)u)v] 
n(n + l)(n + 2)f 3 



-^mm — — Cf (*Sn + 5,^ 



(69) 



where the auxiliary functions read 

6>i(u,<u) = 9(v-u), G 2 (u,v) = 9(-v)9(l-vu), 3 (u,v) = 9(v)9(v -u). (70) 



L n (u,v) 



M?(u,v) 



e(v) fu 



n+l 



V — U \V 

-(— 

m 3 I n + 1 



,n+l 



— U 



n+l 



M 2 >,^) = ^ 



N n {u,v) 



w 3 [n + l 
■ue(w) | 
u(v — u) 1 // 






+ e(v 



n+l' 



u 



f 


1 u n+2 


n + 2 


^n+2 


f 


Mr] 


n + 2 



n+l' 



(71) 



The plus-prescription is defined as 



Gi(u, v)[A n (u, v)]+ = Oi(u, v)A n (u, v) - 6(u - v) I dv'0i(u, v')A n (u, v'). (72) 

It is not difficult to observe that, in multicolour limit, Eq. flBTD is exactly reduced to the equation 
of Ref . |16| , which was the starting point of their analysis. However, we will start from Eq. (|59|) in 
the mixed representation and show in the next section that in the multicolour limit the generalized 
splitting functions can be diagonalized. 
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5 Solution of the evolution equations in multicolour QCD 

This section is devoted to the solution of the evolution equations we derived in the previous 
sections for the twist-3 correlation functions. First of all, we perform an extensive numerical 
study of the exact Eq. fl6"T|). For simplicity we restrict ourselves to the homogeneous case, i.e. 
we discard the quark-mass operator, which is certainly a justified assumption for the light u- 
and d-quark species. The solution we are interested in is given in terms of the eigenvalues and 
eigenfunctions of the anomalous-dimension matrix Z zl l n °f the local operators Z l n calculated in 
Appendix C. The eigenvalue problem we have attacked has no analytical solution; however, the 
diagonalization can be done numerically for moderately large orbital momentum n, e.g. n < 100, 
which is quite sufficient for practical purposes. Second, we provide the analytical solution of 
the generalized DGLAP equation in the multicolour limit of QCD. In this case it reduces to the 
familiar ladder-type equation that holds for the twist-2 operators. It will be shown that such a 
reduction occurs in the limit x — > 1 too. Assembling the two results allows us to construct the 
improved DGLAP-type equations for two-quark twist-3 distributions. Exploiting some examples 
of the evolution for the moments, we confront these two approaches. In particular, we study 
the accuracy of the improved DGLAP equation with respect to the exact evolution for different 
models of gluon light-cone position distribution for the three-particle correlation function at low 
momentum scale. 

5.1 Evolution of the moments 



To obtain the solution of the evolution equation (67) we choose n as a positive integer. In this 
case, as follows from the definition ([58]) of z n (u) the n-th moment is actually given by the following 
linear combination of local operators Z l n (see Eq. ( |G.1| )): 



z" u 



ECtitt^tt'- 1 ^ (73) 



1=1 



so that z n {u) is a polynomial of degree n— 1 in u. Thus the kernel P zz (u, v) possesses n polynomial 
eigenfunctions ef(v): 



dvP2 B (u,v)e?(v) = -\?e?(u), Z = l,...,n, (74) 

where —A" denotes the eigenvalues. These eigenfunctions can be constructed by diagonalization 

„ n 

C k n -l\\dvP n zz {u,v)v n ~ k v k ^ = Y. C n-\zzllu n ~ l u l -\ (75) 



1=1 
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where the anomalous- dimension matrix zzlYk °f the local operators is given by Eq. 
this is a purely algebraic task and we find 



e n M 



Y,C l nAu n - l u l - l El, with {(E n r 1 zzl n E n } M = -X n k S(k-l) 



Actually, 



(76) 



where 5(k — I) is a Kronecker symbol defined by Eq. ( |C.9|) . The spectrum of the eigenvalues Xf 
up to n = 50 is shown in Fig. 0a, together with some examples of the eigenfunctions of the kernel 
P zz for particular orbital momenta. The solution for the moments z n (u) (in the massless case) is 
then expressed in terms of the eigenfunctions and eigenvalues we have found: 



n ( r 

2 n (u,Q 2 )=J2c?(Ql)enu)exp -/ 
i=i ' I Jc < 



3 2 dt a 3 {t) * 

I t 2n l 



(77) 



The coefficients c" (Qq) at the reference momentum squared Ql have to be determined from the 
non-perturbative input z n (u, Qq 



2\. 



lk-l 



t(K)- l (n ^1 ?\ . (1 -'/T- 1 .;" 



cr(Qoy ^~,,- {n _ 1)Uw k-i 



fe=i 



i i ' ^o 

l + w ; \ w=0 



(78 



5.2 Reduction of the evolution equations 

In the large- A^ c limit only the planar diagrams (Fig. 1 a, d) survive and the kernel Pzzi u i v ) 
has two known dual eigenfunctions: 1 and 1 — 2m, so that J duef(u) = 5n + 0(1/N C ) and 
J du (1 — 2u)ef(u) = 5i2 + G(l/N c ), where / = 1, 2 correspond to the lowest two eigenvalues of the 
spectrum shown in Fig. 0a. A straightforward calculation gives the following DGLAP evolution 
kernels: 



Jo 



du 



1 

l-2u 
X-2v 



P zz {x,u,v) = N c 6(x)9{x)< 



+ -x 2 



3 2 

2 



%8{x) 

5{x) 



+ o 



N, 



(79) 



and for the mass-mixing kernels the exact results read 



du 



du< 



1 

l-2u 
\-2v 

1 

l-2u 
l-2v 



P zm (x,u,v) = C F 9(x)9(x 



P zm (x, u, v) = C F 9(x)9(x) 





x(2 — 3x) 



for e, 



for hi 



W) 



As was first observed in Ref. [14| in the context of the chiral-even distribution g2{x), similar 
equations hold true also for the -^--suppressed terms in the x — > 1 limit for flavour non-singlet 
twist-3 evolution kernels. In the present chiral-odd case, we find 



du 



l-2-tt 
1-2-u 



P zz (x,u,v) 



~9(x)9(x) 



+ \8{x) + 0{x c 



+ §S(x) + 0(x 



^.o\ 



+ N C 
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where the N c ■ ■ • symbolize the x — > 1 limit of Eq. (|79| ) . 

The eigenfunctions we have obtained coincide precisely with the coefficients that appear in 
the decomposition of e(x,Q 2 ) and hL(x,Q 2 ) in terms of three-particle correlation functions. To 
observe this explicitly we need relations similar to the ones given by Eqs. ( |i0|) and (|30| ) transformed 
to the mixed representation (158]). Namely, we have 



e(x) 
h L {x) 



—mix) — / duz(x,u), 

x 2 ax Jo 



1 d 



-mix) 



du (1 — 2u)z(x, u), 



x 2 dx Jo 

where we introduce for convenience a new function h L {x), so that hi(x) reads: 

■ :1 d f 1 dy x 2 j 



hrix) 



X 



dy^hiiy) -x 

y 



.,-.Lw/, , , 2 h L (y)- 

ax Jx y y z 



(82) 
(83) 

(84) 



For the massless case the last term on the RHS coincides with the twist-3 part hi,. 

From the observations we have made above, it follows that in the large- N c as well as in the 
large- x limit the twist-3 distributions satisfy the DGLAP evolution equations. By combining 
the large-iV c evolution with the large- x result for the -^--suppressed terms, we can improve the 
accuracy of such an approximation within a factor 5-10, to be compared with the multicolour 
limit taken alone (see Fig. |]b). Thus, the functions e(x,Q 2 ) and hi(x,Q 2 ) obey the following 
improved evolution equations: 




Qs 



dy 

2n Jx y 



Pee(y) e (§) + [i {P mm (y) - Pee(y)} - \Pem{v)£\ ™ (f 

P- h M h L (f ) + [I {Pmm{y) - Pfa{y)} - \P~ hm {y)±] m ( 



*5) 



with 



Uv) = 


--2C F 


V 

y. 


+ ^y+(^-C A )s(y)+0(f/N c ), 


\h(y) = 


--2C F 


V 

y. 


3C 2 Ay+ C C 6 F 4C 3 A ) s ® +0 (y°/ N < 


em(y) = 


= Cfx{2 — x), 


\m{y) : 


= Cpx 


(2 


— 3x). 



id) 



17) 



The evolution kernel for the mass-dependent correlator was already given by Eq. (|63|). In Eq. ( |86| ) 
we have added subleading terms O (y°/N c ) (which we do not specify here) so that the first moment 
of each kernel coincides with the corresponding eigenvalue of the kernel P™ z (u, v). This guarantees 



that the solution for the lowest moments given below, in Eq. (p4|), will be reproduced exactly. Note 
that in the massless case Hl{x) fulfills the same evolution equation as Hl{x). The simplest way to 
verify this is to make the Mellin transform of the corresponding evolution equations. 
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Let us add a few remarks on the momentum space formulation. As we have seen above the 
solution of the evolution equations in the asymptotic regimes is the most straightforward in the 
light-cone position representation. It is by no means trivial to observe the appearance of the 
DGLAP equations in momentum fraction representation. However, we know that the asymptotic 
solution, in coordinate space, is given by the convolution of the three-particle correlation function 
with the same weight function that enters in the decomposition of the two-parton correlators at 
tree level. With this in mind, we are able to check that the integrals 



e(x) 



h L (x) 



d/3'D(x,P') 



'tflifrho-n 



d_ _d_ 



D((3,(3') 



J9) 



taken from Eqs. ([TO) and (f50) neglecting quark-mass as well as twist-2 effects, satisfy the DGLAP 
equations, namely 



e(x) 



a 



-—N c j d(3e((3) { 2 



P 



h L (x) 



a 

Air' 



N c d(3h L {(3){2 



(x-(3) 



(x-f3) 



e° n {x,x-p) 



e° n (x,x-p) 



+ (2-^jQ° 11 (x,x-f3)- l -5((3-x)\, (90) 



+ (2 + 3^j el^x-0) - l -8{(5-x)\ . (91) 



The corresponding anomalous dimensions are 

a 



[e] n 



Air' 



-N r 



mn + 2)-2, E+ 1 - + -±- 2 



^-N c \-2^(n + 2) - 2 lE + \ - -i- 
47T I 2 n + 2 



[h L ]n- 



(92) 
(93) 



Which are exactly the anomalous dimensions 7^ found in Ref. |1J| for e and Hl, respectively, 
with the replacement n — > j — 1 5 . 



5.3 Examples of the evolution 

For the lowest few moments the evolution equation (p7|) can be solved exactly. Taking into account 
the symmetry properties of the quark-gluon correlation function, we obtain the following results 
for the (non-vanishing) first two moments: 

55 [Q 2 dta s (t) 



^ («,«>) = ^(Q S )exp{-!/* 



t 2vr 



5 The difference in the anomalous dimensions is due to an extra power of the momentum fraction x included in 
the definition of the twist-3 correlation functions. 
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(«,«») = ^)e*p{-g/J 



2 f... /y-\ _ ,2(,y2\ J 73 [Q dta s (t) 



t 2vr J 

52 rQ 2 dta s {t) 



where ,z n (Qq) are related to the following matrix elements of the local operators (up to normal- 
ization) 

z l (Ql) = Zl(Ql) = (h\Zl\h) w=Qfj , z 2 (Ql)=Z 1 2 (Ql) = (h\Z 1 2 \h) w=Ql , 
l 2 (Ql) = Zl(Qt) = (h\Zl\h)^ =Ql , z 3 (Ql) = zf(Ql) = {h\Zl\h) w=Ql . (95) 

Thus the Q 2 -dependence for the first moments of e(x) and hi(x) can be predicted uniquely 
since the initial values at the low-momentum scale are given by these moments themselves. For 
larger n, the evolution is sensitive to the shape of the gluon distribution between the quark fields. 
However, as can be traced from Figs. [3] and f|, the dependence on the assumed different toy 
models of the light-cone position distribution at Q% = 1 GeV 2 is quite small: the typical relative 
deviations for the moments of the twist-3 unpolarized and polarized structure functions are 2% 
and 5%, respectively, at Q 2 = 100 GeV 2 . In the calculations we have set the number of flavours 
Nf = 3 and Aqcd = 0.25 GeV. For the "gap" -type (end-point-concentrated) distributions (see 
Figs. [| |] dashed line) the deviation is small with respect to the "coefficient function" -type model 
predictions (1 and 1 — 2m, solid lines), while for the "hump" (end-point-suppressed) distributions 
(dash-dotted line) it is a little bit larger. The accuracy of the multicolour approximation is about 
15-20% at a scale Q 2 = 100 GeV 2 . These numbers are natural and can be expected from the 
discrepancy between the DGLAP anomalous dimensions and the exact lowest two eigenvalues of 
the spectrum (see Fig. |||b). Nevertheless, there is one very important exception from the naive 
expectation. If the initial gluon distribution is strongly suppressed in the end-point region, for 
instance as uK™ -1 )/ 2 ! for u — > 0, then the large- iV c approximation breaks down for large n (polarized 
moments are more sensitive than unpolarized ones). In this case the evolution is not smooth. The 
shape of this function will be turned immediately into the end-point-concentrated one. So, we can 
get rid of this "hump" -type model for a momentum transfer Q ~ 1 GeV and argue that such a 
distribution could not occur in the non-perturbative domain either. It is the most likely to assume 
that the momentum fraction-function Z(x, x') is rather smooth. From the equation 

z n (u) = f dx f dx\xu + x'u) n - x Z{x, x) (96) 

it then follows that z n {u) cannot be strongly suppressed in the end-point region. For instance, if 
Z(x, x') is positive-definite and concentrated in the region < x,x', then z n (u) cannot vanish at 
u = 0,1 unless it is identically zero. 



6 Discussion and conclusion 

In the present paper we have investigated the Q 2 -dependence of the chiral-odd distributions of 
the nucleon e(x) and hi(x). Using the constraint equalities coming from the equation of motion 
and Lorentz invariance, which provide certain sum rules for the structure function, the problem 
is reduced to a study of the renormalization of the multiparton correlators in lowest order of the 
perturbation theory. As a result we construct an exact (taking account of 1/N% effects) one- 
loop evolution in the light-cone fraction as well as in the light-cone position representations. For 
these purposes, we have used two techniques, which employ the light-like gauge for the gluon 
field. Accepting different prescriptions on the spurious pole in the gluon propagator, we were 
able to verify that they do lead to the same results. From the calculational point of view the 
momentum space technique is much easier to treat. However, the coordinate space makes the 
involved symmetries apparent and, as a by-product, diagonalization of evolution kernels is easy 
to handle. We establish the bridge between different formulations of the QCD evolution. It is 
straightforward to obtain the evolution kernels in the light-cone fraction representation, starting 
from the coordinate space and vice versa using the Fourier transform. Using this transformation it 
is straightforward to obtain the coordinate space two-particle evolution kernels for Faddeev-type 
equations with pair-wise particle interaction, which govern the scale dependence for higher twist 
(more than 3) quasi-partonic correlators [ID and to study the diagonalization problem, which is 



more easier to do in the light-cone position representation. We hope to return to this question in 
the future. 

In the multicolour limit as well as for x — > 1 we obtain the ladder-type evolution equations for 
the twist-3 part of the distribution functions. Joining these asymptotics together we construct an 
improved DGLAP equation, which generally has very good accuracy, at the level of few per cent. 
We argue that the observed discrepancy between these reduced equations and the exact evolution 
could occur only for unphysical initial conditions of the latter. To clarify the situation completely, 
it would be helpful to have the low-energy model predictions for the distribution of the gluon field 
in the quark-gluon correlator. 
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A Definition and some properties of 0-functions 



/OO (jry n 

a m H(ax k -l+iOr h . (A.l) 

-oo Mil ,„_-, 



The 0-functions entering the evolution equations in the momentum fraction representation are 
given by the formula 

da 

a'" 

fe=i 
For our practical purposes it is enough to have an explicit form of the functions 

0°(s) = 0, (A.2) 

e°(x) =*(*), (A.3) 

e?^,^) = °(*iM-*2) - °M°(-xi) (A . 4) 

X\ — x 2 
since the other are expressed in their terms by the relations 

e°( x X2 ) = -J??—qO( X2) (A5) 

Xi - x 2 

ei( x X2 ) = -J^e° u (x 1 ,x 2 ) — e° 2 ( Xl ), (A.6) 

Xi - X 2 Xi- x 2 

©22(^1) x 2 ) = -7 1 -^—^Q° u (x 1 , x 2 ), (A.7) 

[xi-x 2 y 

Qi U (x 1 ,x 2 ,x 3 ) = Bn(x2,x 3 ) B?^^!,^), (A.8) 

X\ — X 2 X\ — x 2 
0} n (xi, x 2 , x 3 ) = 6?i(a;2, X3) Q° 11 (x 1 ,x 3 ). (A.9) 

X\ — X 2 X\— x 2 

In the main text we have used the relations 

/ dpp n e° n (p, (3 - x')e° u (x, x-(3) = x' n E n (x, x-x',x) = - \x' n - x n ] O n (x, x-x'). (A.10) 
J n 

PV / 'df3j^pr[e n (t3,p-x) + Q 11 (x,x - /3)] =0. (A.ll) 

B Evolution equations in Abelian gauge theory 

In this appendix we present a pedagogical illustration of the renormalization group mixing problem 
for the redundant basis of correlation functions defined by Eqs. (U)-(|) and Eqs. (|20l) -([25l) for the 
unpolarized and polarized cases, respectively, in the framework of the Abelian gauge theory. Its 
aim is to show the self-consistency of the whole approach we have used as the equations derived 
below satisfy the constraint equalities given by Eqs. flUf), (p6|), (p?|), which are further employed 
to reduce the overcomplete set of correlators to the independent basis of functions. 
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In the calculations of the corresponding evolution kernels, we follow the methods developed in 
Ref. || (see Ref. |21| for a more recent discussion of the RG equations for the time-like twist-3 cut 
vertices). The one-loop Feynman diagrams giving rise to the transition amplitudes of two-particle 
correlation functions into the two- and three-parton ones are shown in Fig. |5]a, b. The last figure 
(c) on this picture is specific of the vertices having non-quasi-partonic form [ID|, that is for e(x) 
and hi(x); it displays the addendum due to the contact term that results from the cancellation of 
the propagator adjacent to the quark-gluon and bare vertices. As an output the vertex acquires 
the three-particle piece. The radiative correction to the three-parton correlators are presented in 
Fig. a, b, c. 

A straightforward calculation yields the evolution equations for the spin-independent case in 
the form 



M(x) 



ex 



a 

2^ 



d(3M((3) { 2 



P oo 



& n {x,x-0) 



_{x-(3) 

^-JdpLp) l^GP u (x,x-0) + i*09-s)J 

P 



+ ^e° 11 (x,x-p) 



(b.i; 



m) {\i*-P) 

0D{p,p)[2 



e° n (x,x-/3) 



+ xOl^x, x - p) + 26?! (a;, x - (3) 



P no 



(x-P) 



e n (x,x-/3) 



+ — -e° ul (x,x-p,x-p + p') 

x — p 



+ 5(P - x) J dp"^Ql n (P^P" - P,p" - P') + 2G° n (x,x - P)\\ (B.2) 



D(x,x') 
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27r' 
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—e(x) - M(x) 
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vv ii IX j X XI 



X 



-e(x') - M(x') 
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Ky iiix.x X) 



+ 
+ 



jn>{ 



, D(x, P') [l3 ', X+ ,f ) e' ln (x' t x' -x,x'-x + If) 
\ [x — x') 



X 



x'-P 



7 [D(z -x' + p', p') - D(x, x')]e° n (x', x' - pi) 



I dp(D(P,x') ^, f + . x) Q° ni (x,x - x',x- x' + (3) 



+ 



(x' — x) 
-^—[D(P, x'-x + P)- D(x, x')]<d° u (x, x-P))- -D(x, x') 



(B.3) 



where the dot denotes the derivative with respect to the UV cutoff ' = A 2 <9/<9A 2 and the plus- 
prescription is defined by the equation 



P 



{x-P) 



e° u (x,x-p) 



jJL^& n {x, x-P)-6(P-x)Jdp"—^e° n (P", P"-P). (B.4) 
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In the same way we can immediately obtain the set of evolution equations for the spin- 
dependent distributions (p0|) -(p5|): 

P + x 



M(x) 

h(x) 



rJ d ^{ 2 



e° n (x,x-f3) 
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2^ 



d(3h 1 {(5){2 
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(x-{3) 
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-e° n (x, x -p) 



(B.5) 
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The anomalous dimensions calculated from the evolution equation for the distribution h\(x) co- 
incide (up to the colour group factor Cp) with the result of Ref. p2 |. 
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By exploiting the relation provided by the equation of motion and Lorentz invariance, we can 
easily verify that the RG equations thus constructed are indeed correct and the renormalization 
program can be reduced to the study of logarithmic divergences of the three-parton (Z(x,x'), 
Z(x,x')) and quark mass (M(x), M(x)) correlators in perturbation theory. 

C Local anomalous dimensions 

In this appendix we pass from the evolution equations for correlators to the equations for their 
moments and, in this way, the anomalous-dimension matrix for local twist-3 operators. 
We define the moments as follows: 

F n = Jdxx^x) for any two-particle correlator, 

Z l n = J dxdx'x n - l x a - l Z(x,x'). (C.l) 

In the language of operator product expansion these equalities specify the expansion of non-local 
string operators in towers of local ones, namely: 

al— 1 an— I 

= ^(OKiD+J'-'gG+^O)^ f ' J («? + )-V(o), 



d n . ,, , s , m 

75 

The inverse transformations to the non-local representation are given by 

. m \n oo \ 

777 I 77 I ' 

n=0,m=0 " L - IL ' n=0 



M n = t"—M(\)\ x=0 = -mi + \ I (i£>+rV>(0). (C.2) 



Z{\^)= E Hr^i-^r^-jZ^+i, M(X) = y £(-i) n ^rM n . (C.3) 

ml nl ~k n\ 



Now it is a simple task to derive the algebraic equations for the mixing of local operators under 
the change of the renormalization scale from the evolution equations (|35|)-(|39D. They are 



Mn = ^MMl n M n , (C.4) 



• l a 

Z1T 



n n n - l+1 ± zmH\ M n + £ zzl iz k n , (C.5) 

fe=l ) 



where the anomalous dimensions are given by the expressions 

mm7™ = ~Cf (S n + S n+2 ) , (C.6) 

zm7 " = /(/ + l)(/ + 2)' (C ' 7) 
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Here we have used the following step functions 
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Z -> n- l + \ 
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1, i > J 
0, t < j 
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m\{n—m)\ ' 



(C.9) 



The plus and 



as well as the convention S n = YX=i % anc ^ ^ ne binomial coefficients C™ 

minus signs in this equation correspond to the functions e and hi,, respectively. These analytical 



expressions coincide with the result of Ref. [15 
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Figure 1: The one- loop renormalization of the three-parton correlation functions. 
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Figure 2: The spectrum of the eigenvalues A" for the evolution kernel P™ z defined in (JB9|) is 
shown in (a). In (b) the relative deviation 1 — A"/(— 7") (in %) for the lowest two eigenvalues 
of the spectrum, i.e. I = 1,2, is plotted: the solid (dashed) line 7* (7^) corresponds to the n-th 
moments of the improved kernel P ee (P^) defined by Eq. (|56D; the dash-dotted (dotted) line is the 
multicolour approximation for P ee (P^). In the improved approximations, subleading terms were 
taken into account to reproduce the first two eigenvalues exactly. Eigenfunctions of the kernel 
Pzz are shown for n = 4 in (c) and for n = 30 with I = 2, 3, 6, 29 in (d). 
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Figure 3: The predictions of the improved evolution equation for the moments [e] n , normalized 
to 1 at Qq — 1 GeV 2 , are shown in (a) for n = 1 (solid line), n = 2 (dashed line), n — 10 
(dash-dotted line) and n = 100 (dotted line). In (b) three different models for the gluon light- 
cone position distributions at Ql = 1 GeV 2 for n = 10 are presented; they are evolved up to 
the scale Q 2 = 100 GeV 2 in (c). The relative deviation 1 - [e(Q 2 )]° x /[e(Q 2 )]jf 1 (in %) from the 
improved DGLAP equation of the exact evolution for the assumed gluon distributions is shown 
in (d). The dotted line corresponds to the relative deviation 1 — [e(<5 2 )]^ c /[ e (Q 2 )]n° ef with respect 
to the multicolour approximation, where "coef ' refers to the gluon distribution that is equivalent 
to the corresponding coefficient function. 
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Figure 4: The predictions of the improved evolution equation for the moments [hi] n , normalized 
to 1 at Qq = lGeV 2 , are shown in (a). Here the solid line and dashed line represent n = 2 and 
n = 3, respectively. The further description is the same as in Fig. [j], except that n = 20 in (b)-(d). 
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Figure 5: One- loop radiative corrections for two-particle correlators in the Abelian gauge the- 
ory. The fermion propagator crossed with a bar on diagram (c) shows the contraction of the 
corresponding line into the point. 
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